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ABSTRACT 



We study the possible magnitudes of CP and lepton-number-violating quantities in specific 
GUT models of massive neutrinos with different Abelian flavour groups, taking into account 
experimental constraints and requiring successful leptogenesis. We discuss SU(5) and flipped 
SU(5) models that are consistent with the present data on neutrino mixing and upper limits 
on the violations of charged-lepton flavours and explore their predictions for the CP-violating 
oscillation and Majorana phases. In particular, we discuss string-derived flipped SU(5) models 
with selection rules that modify the GUT structure and provide additional constraints on the 
operators, which are able to account for the magnitudes of some of the coefficients that are 
- ' often set as arbitrary parameters in generic Abelian models. 
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1 Introduction 



In recent years, several experiments have provided convincing evidence for neutrino oscillations 
PQ [2j [31 HI |5] , together with important information on the neutrino mass differences and mixing 
angles [6]. It has now been established that the atmospheric neutrino mixing angle, #23, is 
large, with a central value ~ 7r/4, so that 0.31 < sin 2 #23 < 0.72. The solar MSW angle #i 2 [Z] 
is also large, but not maximal, with sin 2 # 12 ~ 0.3. Upper limits from the CHOOZ experiment, 
in particular, establish that the third mixing angle, #13, must be small, with a present upper 
limit of sin 2 # 13 < 0.1 [8]. Finally, experiments have established the following ranges for the 
mass differences: Am 2 4m ~ 2.5 10~ 3 eV 2 and Am 2 oi ~ 8 10~ 5 eV 2 [6]. 

These important results leave several questions open for model builders to consider: 

- Is the atmospheric neutrino mixing indeed maximal? 

- How large is the solar neutrino mixing, and can it be related to mixing among quarks? 

- What is the pattern of masses for the neutrinos? 

- What is the magnitude of 9 13 7 

- Is there significant CP violation in the neutrino and charged-lepton sectors? 

In this paper we study, in particular, the answers to the last two questions that are provided 
by models that propose different answers to the first three questions. 

The existence of neutrino masses and oscillations may have further consequences. For example, 
the decays of heavy neutrinos may lead to leptogenesis, which would impose additional con- 
straints on CP- violating phases and on the magnitudes of certain neutrino Yukawa couplings [9]. 
These should be large enough to generate a sufficient net lepton number, but not above the 
value that would erase any generated asymmetry through equilibriation effects. Also, mixing in 
the charged-lepton sector is enhanced by radiative corrections in supersymmetric models that 
mix the sleptons, contributing via loop diagrams to rare decays and flavour conversions [TmfTT]. 
The rates for such processes mediated by sparticles may be large and close to the current ex- 
perimental bounds in many supersymmetric models. Moreover, the possibility of CP-violation 
is crucial in this analysis [12] and additional CP-violating quantities may be observable in the 
charged-lepton sector. The model parameters that would be responsible for leptogenesis, flavour 
and CP violation in charged-lepton processes are different from those measurable in neutrino 
oscillation experiments, and also motivate this paper. 

As we discuss below, the predictions for these additional parameters vary in different models. 
However, there are also several correlations. A combined analysis, using all available processes, 
could give significant constraints on models, and even exclude certain classes of models. The 
joint study of phases, CP and lepton-flavour violation is a central issue in this programme, and 
its better understanding is the central goal of this work. 

The fact that the fermion mass matrices exhibit a hierarchical structure suggests that they are 
generated by an underlying flavour symmetry for which the simplest possibility would be an 
Abelian family symmetry. In such a case, one can parametrize the charges of the Standard 
Model fields under the symmetry as in Table 1. The Higgs charges are chosen so that the terms 
fzf^H (where / stands for a generic fermion, the subscript 3 denotes the third family, and 
H denotes either Hi or H 2 ) have zero net charge whereas other terms have non-zero charges. 
Then, in general only the (3,3) element of the associated mass matrix will be non-zero as long 
as the U(l) symmetry remains unbroken. The remaining matrix elements may be generated 
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when the U(l) symmetry is spontaneously broken [T3J [H] by fields 9, 9 that are singlets of 
the Standard Model gauge group, with U(l) charges that are in most cases taken to be ±1, 
respectively Here we assume them to have equal vacuum expectation values (vev's). The 
suppression factor for each entry depends on the family charge: the higher the net U(l) charge 
of a term fifjH, the higher the power n in the non-renormalisable term fifjH (jj-^j or (|j) n ] 
that has zero charge, where M is a high mass scale that is assumed to be much larger than 
(9) = (9). The f 3 f 3 H couplings are not suppressed, because they have no net charge. 



U(l) 



Qi u\ d\ U e\ 



Ho 



a; 



Pi 7i 



di 



-a 3 - p 3 -a 3 - 73 



Table 1: Notation for possible U(l) charges of the various Standard Model fields, where i is a 
generation index. 



Whilst such a family symmetry provides a promising origin for a hierarchical pattern of fermion 
masses, in order to go further it is necessary to specify the charges of the quarks, charged leptons 
and neutrinos [15] . In principle, one can fit the observed masses and mixing angles by choices of 
flavour charges that are not constrained by additional symmetries. However, the structure of the 
Standard Model is suggestive of an underlying unification relating quark and lepton multiplets, 
as occurs in various GUT models. In such a case, it is important to try to understand the 
different mass and mixing patterns and other correlations predicted by different GUTs, and the 
resulting patterns of CP and lepton-flavour violation in particular. This is the primary focus 
of this paper. Since particles in the same multiplet of a GUT group have the same charge, the 
predictions and correlations between observables will in general be different for different GUTs 
and charge assignments. 

A crucial question when studying flavour symmetries is whether these are indeed Abelian, as 
described above, or non- Abelian. In the case of an Abelian flavour symmetry, taking predic- 
tions for the masses and mixings only from the flavour and GUT structure and assuming that 
unknown numerical coefficients are generically 0(1) leads in the simplest and apparently more 
natural neutrino models to large mass hierarchies and unacceptably small solar mixing, which 
is correlated with the hierarchies of charged- lepton masses |16j . Large solar mixing, as required 
by the experimental data, has therefore to arise either from the right-handed neutrino sector 
(as may arise, for instance, in see-saw models with dominance by a single right-handed singlet 
neutrino, and other zero-determinant solutions), or by imposing more U(l) symmetries and/or 
introducing more fields. However, this strategy introduces additional model dependence and 
loses predictivity. On the other hand, large atmospheric mixing is naturally predicted in a wide 
class of Abelian models, and non- Abelian flavour symmetries typically predict naturally large 
angles for both solar and atmospheric neutrinos [TT1 [TS] . 

Since there is no charge quantisation in Abelian groups, contrary to the non-Abelian case, 
the symmetry cannot determine alone the numerical values of the mass-matrix elements and 
hence the mixing angles and phases. Nevertheless, information on the phases and values of the 
coefficients can be inferred phenomenologically, and individual Abelian models may be used to 
make interesting qualitative predictions, just as they have for the mixing angles, based on the 
powers of the small parameters (9) /M, (9) /M that are allowed by the symmetries. However, any 
predictivity requires the existence of stable solutions where a small change in the parameters 
does not change drastically the predictions. Non-Abelian symmetries are more predictive in 
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this respect. However, once the symmetry is broken (as is required to generate large lepton 
mass hierarchies) this predictivity is to a certain extent lost. Even in this case, therefore, the 
phases are derived by a combination of theoretical and phenomenological considerations, again 
in a mo del- dependent way. 

The fact that neutrinos have masses and mix with each other like quarks also implies that, in 
principle, we can expect non-negligible CP- violation in the neutrino sector. This may manifest 
itself in several different ways in low-energy neutrino physics, including the Maki-Nakagawa- 
Sakata (MNS) oscillation phase 5 and two possible Majorana phases. Additional phases arise in 
extensions of the light-neutrino sector to include either heavy singlet neutrinos and / or charged 
leptons. For example, even the minimal three-generation see-saw model has 6 independent CP- 
violating phases, and leptogenesis is independent of the light-neutrino phases, in particular. 
On the other hand, some of the phases on which leptogenesis does depend may in principle be 
observable in the charged-lepton sector, in the presence of low-energy supersymmetry. As we 
discuss below, these observables may take very different values in classes of models that fit the 
present neutrino data equally well. 

The structure of this paper is as follows. In Section 2 we discuss the neutrino and charged- 
lepton observables for which we investigate the predictions of various GUTs with Abelian 
flavour symmetries. The specific GUT models and their predictions are discussed in Section 3. 
Our conclusions are set out in Section 4, and appendices present details of the specific GUTs 
studied. 



2 Observables in the Lepton Sector 
2.1 Neutrino Observables 

The mixing matrices K,r>, Ur ;U that yield physical CP and flavour violating observables are 
those making the following diagonalisations: 



where Y e ^ stand for the Yukawa couplings of neutrinos and charged leptons respectively, and 
Mrr is the heavy Majorana neutrino mass matrix, which we assume to be three-dimensional 
(for a review see [H]). The effective neutrino mass matrix is: 




(1) 
(2) 
(3) 
(4) 



1 




(5) 



m e ff m D 



M RR 



In terms of the above matrices, the Maki-Nakagawa-Sakata (MNS) matrix becomes 



u MNS = u = v e ul 



(6) 
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and can be parametrized as: 



U = diag(e i5 % e*", e^.V.diagie-^ 2 , e" i<fe/2 , 1) 



where 



V 



C12C13 



■Sl2 c 13 



-C23S12 - S 2 3Sl 3 Ci2e 1 ' 5 C23C12 - S 23 S 13 Si 2 e l5 



Si3e 
S23C13 



;<5 



\ S23-S12 - C 2 3Si3C 12 e l5 -S23C12 - C23Si 3 Si 2 e i ' 5 C23C13 



(7) 



(8) 



and Cij and s^- stand for cos 8 ij and sinfly, respectively. In this formalism, the three neutrino 
mixing angles and the six CP-violating phases are given by 



#13 = arcsin (If/isl), #12 = arctan 



\U- 



12 



?23 



arctan 



1 17: 



23 



33 I 



I J n + Ci 2 C2 3 Si 3 



5 M = arg(C/ 2 3), $t = arg(C/ 33 ), 5 = -arg 
<5 e = arg(e i5 f/i 3 ), 0i = 2arg(e i5e f/i*i), 2 = 2axg( e ^£/* 2 ) . 



V 



S12S23 



where i,j = 1, 2, 3 and i ^ j. A useful measure of the amount of CP-violation in the oscillations 
of light neutrinos is the Jarlskog invariant 

Jcp = ^\im(u* 1 u 12 u 21 u; 2 )\ = ^\im(u^ 1 u 13 u 31 u; 3 )\ 

I I 

= -|/m(C/ 2 * 2 C/ 2 3^32^3*3)l = 2l c i2 c i3 c 23sin5s 12 si 3 s 2 3| ■ (9) 



It is clear that different models that reproduce the experimental values of #i 2i23 and the differ- 
ences in light neutrino masses-squared may well predict different amounts of CP-violation in 
neutrino oscillations ([9} and elsewhere. So far, experiment has provided only an upper bound 
on the possible magnitude of #13, and no information on 5. The magnitude of CP- violation in 
neutrino oscillations may even be zero, for example in models with texture zeroes in the (1,3) 
entries. As we see later, specific models may connect #i 3 and 5, and also provide connections 
to other observables. We aim here at the differentiation of possible flavour models according 
to their predictions for observables in the neutrino sector - what are the expected magnitudes 
of 9i3,8 and 02,3? - and elsewhere 0. We now discuss some of the other observables in more 
detail. 



2.2 Leptogenesis 

The textures of neutrino Yukawa couplings may be constrained by the requirement of successful 
leptogenesis [2T], which requires obtaining the correct magnitude of lepton-flavour violation 

order to calculate mixing angles and CP- violating phases in an automated way, we analyze the different 
models using the package REAP [2U] , 
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(LFV) in the decays of heavy, right-handed Majorana neutrinos via a difference between the 
branching ratios (BR) for heavy-neutrino decays into leptons and antileptons: 



BR{Nl -> $ + £) + BR(N C L -» $ + £). 

Since non-perturbative lepton- and baryon-number- violating interactions mediated by sphalerons 
are in thermal equilibrium up to the time of the electroweak phase transition, a non-zero lepton 
number gives rise to a non-zero baryon number, by sharing the lepton asymmetry AL ^ with 
a baryon asymmetry AB ^ 0. 

In order to avoid washout of the initial decay asymmetry by perturbative decay, inverse decay 
and scattering interactions, the model must satisfy an out-of-equilibrium condition, namely 
that the heavy-neutrino decay rate is smaller than the Hubble parameter H at temperatures of 
the order of the right handed neutrino masses. The tree-level width of the heavy neutrino Ni 
with mass Mi is: T = [(A^A) aj '871} Mi, which should be compared with the Hubble expansion 
rate H « 1.7 gl /2 T 2 /M P , (where gf SSM « 228.75, gf M = 106.75 and M P is the Planck 
mass), leading to the requirement 

(AU)ii M P < M^ 



UTigl' 2 



Here, A = Y' v denotes the neutrino Yukawa matrix in the basis where the Majorana masses Mi 
are diagonal. This condition may be implemented more accurately by looking in detail at the 
Boltzmann equations, but this formula is sufficient for our purposes. 

The CP-violating decay asymmetry €j in the decay of a heavy-neutrino flavour j arises from 
the interference between tree-level and one-loop amplitudes, and is 





e,- = ;t Im 

where the kinematic function 

f(y) = Vv 

The first term in f(y) arises from self-energy corrections, while the second and third terms 
arise from the one-loop vertex. For models with degenerate Majorana masses, one expects a 
resonant enhancement of the lepton asymmetry, since in this case y ~ 1 and 

m ~-2(\M^-\ Mt \y 

This will be manifest in some of the examples below. 

The lepton asymmetry produced by the decays of the heavy neutrinos is in general diluted 
partially by lepton-number- violating processes [21]. The washout factor Ki is approximately 



, . /lO" 3 eV\ / fhi \~ os , , 

Kj(fti) ^ . 3 l 0g _L_ , (12) 
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where 

rh^iM^M^u/mi (13) 
and the primes denote the basis where the right-handed neutrino mass matrix is diagonal. 
Then, for a supersymmetric model, the lepton and baryon asymmetries are given by 

UL > ~ 0.8 x 1(T 2 Ki €i (14) 



" J ' M ~ -2.8 x i (15) 



3 / A' 

'-) ~ -2.8 x lO" 3 ^ x . 
s ) A 

•na ^ -0.02/^e, -, (16) 

where the entropy of the co-moving volume s is given by s = (2/A5)g^n 2 T 3 , and g* = 228.75 
refers to the effective number of relativistic degrees of freedom contributing to the entropy in 
our supersymmetric model. 

In the above relations, A is a dilution factor that appears due to entropy production from 
symmetry breaking and an inflationary phase. This effect has been studied in [22j, in the 
context of the breaking of SU{5) x U(l) when a singlet field $ gets a vev. In this case, the 

dilution factor is [23] A = 4f^4(f^) 3 v m 3 ] 2 i -. Here, the $ decay rate is given by 



s(Rd v ) V Rdr, , ,„ s . ( , , 



T$ = a$ tjMfz , V is the scale where the vev of the Higgs 10 and 10 break the flipped- SU (5) 
group, m n is the inflaton mass, and msusY is the supersymmetry breaking scale. 



Successful leptogenesis would require 

T] B = (6.15 ±0.25) x 1(T 10 . (17) 

Clearly, Ki/A < 1 in general and a model with too small a value of would not be valid. 
However, as we discuss below, this is not the case in the models we study, which generally 
require A ^> 1 in order to provide a successful scheme of leptogenesis. An alternative to 
invoking a large dilution factor would be to adopt a wider range of neutrino Yukawa couplings 
and Majorana masses that can lead to consistent solutions. 

The impact of flavour effects in leptogenesis has also been considered in |26j and more recently 
in [27], and may affect the results by factors of ~ 2. Since the conclusions of our paper are 
not affected by the possible presence of such factors, we will not proceed with more detailed 
considerations in this direction. 



2.3 Charged-Lepton-Flavour Violation 

In the context of low-energy supersymmetry, mixing in the neutrino sector also generates mixing 
in the sleptons via loop corrections, contributing to rare decays and lepton-flavour conversions. 
We evaluate these effects in the context of the CMSSM, where the soft supersymmetry-breaking 
masses of the charged and neutral sleptons are assumed to be universal at the GUT scale, and 
may be written in diagonal form with a common value mo- Off-diagonal entries in the slepton 
mass matrix m| are then generated radiatively by the renormalization-group evolution from 
the GUT scale Mq UT . 
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The branching ratios for LFV decays can be described well by a single-mass-insertion approxi- 
mation [281129] : 

BR{U^l n ) « ^(m ,Mi /2l /i)|m| J 2 tan 2 /?, (18) 

where is a function of the soft supersymmetry-breaking masses fixed at the high scale 
Mgut, and % 7^ j are generation indices. In the commonly-used approximation for solving 
the renormalisation-group equations of a single intermediate right-handed neutrino threshold, 
the mass corrections are related to the neutrino Yukawa couplings by [3D] 

2 _ <? lk I A J. 1 A (I \ l<? 3 k \* 



rrif = k 



^Y„(At + A£ k )(Yj ")*, (19) 



k 



where Y v is the Dirac neutrino Yukawa matrix in a basis where both the heavy-Majorana and 
charged-lepton couplings become diagonal, and 

K = -6m 2 - 2A 2 , At = ln(M GC7T /M 3 )/167r 2 , A£ k = ln(M 3 /M fc )/167r 2 . (20) 

where Mj with i = 1,2,3 denote the Majorana neutrino masses assuming Mi < M 2 < M 3 . 

This result is not very accurate, but it is a useful approximation for obtaining analytical ex- 
pressions for lepton-flavour- violating decay rates [31]. In the analysis that follows we assume 
for simplicity that A Q = 0. If A 7^ 0, this parameter would also be similarly renormalized, via 

« ~i^ AoYei{V ^ )ijlog ^- (21) 

The extra suppression due to Y e and the form of k indicate that, unless Aq becomes much larger 
than mo, our conclusions will be qualitatively unchanged. 

Another lepton-number- violating observable is /i — > e conversion on a nucleus, which has a rate 

w 5.6 x 10~ 3 BR(fi -»■ erf). (22) 

This process is very interesting, despite the relative suppression by about two orders of mag- 
nitude, because of the accuracy possible in future measurements of this process. Within the 
framework discussed here, a similar suppression is expected for the decay /x — > 3e, 

r(/i + — ► e + e + e _ ) , 

— k6x 10~ 3 . 23 

r(/i+^e+ 7 ) 1 ^ 

However, the present experimental bound on this branching ratio is relatively weak, and the 
prospects for significant improvement are more distant. It should be noted, though, that this 
decay does offer the possibility of observing a T- violating asymmetry, which might be observable 
if \i — > e7 occurs with a rate close to the present experimental upper limit. The flavour-violating 
decays r — » ey and 7x7 are also potentially interesting. They are governed by formulae similar 
to those above, and later we present some results for them. 
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3 GUTs with Abelian Flavour Symmetries 



As mentioned in the Introduction, we focus in this work on Abelian flavour symmetries, since 
they are simple and arise naturally in a wide class of models. Moreover, understanding Abelian 
symmetries is a first step towards the understanding of non-Abelian symmetries. Indeed, by 
combining two or more Abelian symmetries and introducing more than one field whose vev 
determines the expansion parameter of the mass matrices, one can simulate to some extent the 
picture that one would obtain from a non-Abelian structure. 

3.1 SU(5) 

The first GUT group we analyze is SU(5), whose minimal matter field content is three families 
of (Q,u c ,e c )i G 10 representations of SU(5), three families of (L,d c )i e 5 representations, and 
heavy right-handed neutrinos iV£ = vr in singlet representations. Only two heavy iV£ = Vr 
fields are needed to provide two non-zero masses for light neutrinos, as required by experiment, 
and more than three are present in some models. However, in what follows, we assume that 
there are also just three heavy iV£ = vr fields. This model has the following properties that 
are important for our analysis: 

(i) the up-quark mass matrix is symmetric, 

(ii) the charged-lepton mass matrix is the transpose of the down-quark mass matrix, which 
relates the mixing of the left-handed leptons to that of the right-handed down-type quarks. 
Since the CKM mixing in the quark sector is due to a mismatch between the mixing of the 
left-handed up- and down-type quarks, it is independent of mixing in the lepton sector. In 
particular, in SU(5) the large mixing angle that is observed in atmospheric neutrino oscillations 
can easily be consistent with the observed small Vckm mixing. Most importantly, 

(iii) it is clear from the SU(5) representation structure that the Abelian flavour charges of the 
fermions (Q, u c , e c )j in the same 10 representation must be identical, as must the charges of the 
(L, d c )i in the same 5 representation. This is the type of correlation that we seek to test by 
examining model predictions for flavour and CP violation in the generalised lepton sector. 

There are several possibilities for the Abelian flavour charges, that are motivated by theoretical 
considerations such as anomaly cancellation [32] as well as phenomenological arguments. Within 
this framework, one may address questions such as: 

- How close to maximal is the atmospheric mixing? 

- How large is the solar mixing? 

- How large are 6*13 and 5? 

- At which level is mixing controlled by the hierarchies of charged-lepton masses, and how 
strong is the influence from the heavy right-handed neutrino sector? 

- What statements can be made about other observables? 

There are several viable sets of textures that we analyze for answers to these questions. Tak- 
ing into account the above-mentioned constraints arising from the SU(5) multiplet structure 
(symmetric up-quark mass matrix, and Me the transpose of Md OW n), the mass matrices are 
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constrained to the forms 



M u oc 



Jx+b\ 



where a 



e |x+6| 




e |26| 


e \b\ 


e \b\ 


1 


Qlb 


= Q 



M^d own OC 



l\y+*\ 



g|a+6| 



#1 



Mo oc 



g|x+a| 
g|x| 



gl»l 

g|a+fe| g|a| 
#1 1 
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.r 



(24) 

Qi — Q|- We note that, by assumption, 



we have assumed that the flavour charges of the two supersymmetric Higgs fields are the same. 
Since we know that the (3,3) entry in the charged-lepton mass matrix is the largest one, we 
have simplified our considerations by taking a zero Higgs flavour charge 0. 



3.2 An example with large tan/3 

An interesting possibility is that maximal atmospheric mixing arises from the charged-lepton 
sector. In this case, the (2,3) and (3,3) lepton entries are comparable, and a = 0. The resulting 
mass matrices have been discussed in the literature, e.g., in |16j . The choice \b\ = 2 leads to the 
correct m s /mb ratio, while requiring correct (1,2) quark mixing fixes \x\ = 3 in the down mass 
matrix and \y\ is fixed by the down and charged-lepton mass hierarchies. Thus, one obtains 
finally mass matrices of the form [16, 33J: 



e 6 e 5 e 3 



M u oc e 5 e 4 e 2 , M down oc 



e 3 e 2 



( ^ 4 


e 3 






/ 


e 4 


e 3 




e 3 


e 2 


?) 


, M^oc 




e 3 


e 2 


!) 


V e 


1 






I 


e 3 


e 2 





(25) 



where a single expansion parameter has been used to reproduce the fermion mass hierarchies 0. 

As is clear from the form of the charged-lepton mass matrix, unlike what happens for the 
atmospheric mixing, the solar mixing cannot be obtained from the charged-lepton sector. In this 
case, it has to arise from the neutrino sector, exploiting the see-saw structure. It is interesting 
to observe the a priori large value of #13 given by the charged-lepton sector - indeed, #13 may 
even be too large - but the overall (1,3) lepton mixing can be reduced either by a numerical 
coefficient, or by a cancellation between the charged-lepton and neutrino mass matrices. In fact, 
comparing the textures for maximal and non- maximal mixing, we see that, due to the different 
magnitudes of the neutrino mixing angles entering in the Jarlskog invariant, one can expect 
in principle different inter-correlations and predictions for CP-violation, even before looking in 
detail at the heavy right-handed neutrino sector and its influence on m e ff. 

In general, of course, both the mass structure and the mixings of neutrinos are more compli- 
cated, because of the heavy Majorana masses of the right-handed components. We assume that 
these arise from terms of the form z/^i/^E, where £ is an SU(3)x SU(2) xU(l)-invariant Higgs 



2 We look below also at examples with small tan (3. 

3 The conditions for anomaly cancellation give further insight into the possible textures, but detailed model 
building goes beyond the scope of this paper. 

4 We cannot obtain the (1,2) quark mixing from the up sector, as this would lead to an unacceptably large 
mass for the up quark. 

5 This example gives an up-quark mass that tends to be somewhat too high, but this defect can be remedied 
in specific models. 
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field with Iyy = and a non-zero flavour U(l) charge 0. The possible choices for the £ charge 
give a discrete spectrum of forms for the Majorana mass matrix M VR . In the case of a single £ 
field and assuming a zero Higgs charge and left-handed lepton charges as above, we have: 



m 



D 



oc 



?\y+ni\ e \y+n 2 \ e \y+n 3 \ 

]a+n 2 \ Aa+n 3 \ 



\a+m 



-\ n l\ 



|n 2 | 



M RR oc 



z\2ni+a\ p\ni+n2+<r\ 



:\ni+n2+cr\ 



\2n 2 +a\ 



g|ni+"3+o-| g|n2+n3+cr| 



g|ni+"3+o"| 
g|n2+rt3+o-| 
g|2n 3 +(r| 



(26) 



(27) 



where the charges n« are the U(l) charges of the right-handed neutrinos, and a is the U(l) 
charge of the field S. Unlike what happens for the charged leptons, where the (3,3) entry is 
the largest one, the large entry in Mrr can be in any position in the matrix, depending on the 
relative charges of E and the right-handed neutrinos. 

We recall that the effective light-neutrino mass matrix is given by 

„ 1 



rn 



m 



D 



M. 



-m 



RR 



v T 
D ■ 



Its diagonalization relative to the charged-lepton mass matrix (which is the transpose of the 
down-type quark mass matrix) determines the MNS mixing matrix. In the case of a zero flavour 
charge for the Higgs, the mass matrix is determined by the flavour charges of the left-handed 
neutrinos, which in SU(5) are the same as those of the right-handed down quarks. The effective 
light-neutrino mass matrix then has the structure 



rn 



eff 



■oc 



e\y\ 



1 



z \2a\ 



(2f 





( 


e 2 


6 


P 




/ 




g|l±na| 


g|l±n3 


m eff oc 




I 


1 




, m u D oc 






g|n 2 | 


gl"3| 




K 


I 


1 






I 




^\n 2 \ 





From the values of a and y obtained from the charged-fermion mass hierarchies, we then con- 
clude that 

/ e 2 e e \ ( e |1±ni1 e |1±n21 e |1±ri31 \ 

(29) 

which can potentially lead to large solar and atmospheric neutrino mixing (the expansion 
parameter of the Dirac neutrino mass matrix is in principle similar to that for the up quarks, 
since these particles couple to the same Higgs field, but there can be deviations from this). 

A working example along these lines can be obtained by the following choice of charges: 

rii — 2, n 2 — — 1, n 3 = 1, a — -1. (30) 

Using the indicative choices of the coefficients a?- shown in Table [21 we have calculated the 
relevant observables shown in Table [3j The expansions in e are obtained by using similar tech- 
niques as in Ref. [31]. Altering the coefficients a?- would affect the numerical values appearing 
in Table El but not the powers of e. 



6 Note that E is a singlet that enters only in the heavy Majorana mass textures. It is not the same as the 9 
field that generates the light fermion masses. 
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Parameters in an SU(5) model with large tan (3 


Charged leptons 


a\ 2 = 0.6, a\ 3 = 0.9, a\ 2 = 1.2, a e 23 = -0.5e^ 3 , a e 31 = 0.7, a% 2 = 0.6, a% 3 = 0.4 


m u D 


a\ 2 = 1.3, = -1.3, a\ 2 = 0.7, a£ 3 = 1.8e^ 5 , a| 2 = 0.7, a% 3 = 0.5 


M RR 


a& = l,a£ = 1.8 



Table 2: Choice of coefficients that reproduce the fermion data for an SU(5) model with large 
tan/3. Coefficients not listed in the Table are set to unity. 



Observables 


Series Expansions 


Numerical Values 


m U3 (eV) 


0.23e + 0.35e 3 + 0.66e 5 


0.05(0.047) 


m Ul / m V3 


0.85e-0.51e 2 -1.92e 3 


0.13(0.17) 


m V2 /m V3 


0.84e + 0.52e 2 - 1.93e 3 


0.17(0.17) 


Mt/Ms 


1.8e - 1.21e 2 - 0.34e 3 


0.31(0.36) 


M 2 /M 3 


1 + 1.21e 2 + 3.62e 3 


0.95(1) 


m Uee (eV) 


O.obe — 2.9oe + lz.94e 


0.0043(0.0044) 


a 

^23 


u.yz — i.oye — i.ooe 


u.oo^u.yz ) 


U1 9 
1Z 


0.78 - 1.43e + 3.62e 3 


0.53(0.78) 


#13 


0.96e-2.92e 3 + 3e 5 


0.17(0.19) 


,5 


3.72-0.18e + 6.1e 2 


4.3(3.72) 


J op 


0.093e-0.26e 3 -0.57e 5 


0.017(0.019) 


01 


5.6-3.43e+19.3e 2 


5.71(5.6) 


02 


2.46 + 2.68e + 19.3e 2 


3.93(2.46) 


ei 


0.086e 3 - 0.04e 4 + 0.3e 5 


4.1 x 10~ 4 (0.69 x 10~ 4 ) 


£2 


0.005e 3 + 0.02e 4 + % 02 £™f 


0.007 (0.0048) 


£3 


0.005e 3 0.02e 4 + °f^f 


0.0063 (0.0048) 


Vb x A 




3.16 x 10^(2.46 x 10" 7 ) 


16*" x m 2 


13. le- 1.42e 2 - 60e 3 


2.2(2.6) 


£?i?(/x -> e 7 ) x (ta ^ )2 




5.3 x 10^ n (7.4 x 10~ u ) 


x m? 


11.5e + 0.84e 2 - 41e 3 


2.2(2.3) 


BR(r -> e 7 ) x ( ^ )2 




4.4 x 10^ n (5.7 x 10~ n ) 


167r ' x m? 


12. 7e 2 - 2.2e 3 - 51e 4 


0.43(0.51) 


Bf2(r -> /i 7 ) x (tan/3)2 




2 x 10~ 12 (2.8 x 10~ 12 ) 



Table 3: Values of observables predicted in the SU(5) model with large tan/3. The rate of 
convergence of each expansion in e can be judged from the relative magnitudes of the expansion 
coefficients and by comparing the exact numerical value of the observables with that obtained by 
keeping only the dominant term in each expansion (given in parenthesis except for e 2 and e 3 , 
where the number in parenthesis corresponds to the full fraction). 



In this example, Y^Mqut) ~ 0.6, which can be compatible with tan/3 from 40 to 50 (taking 
into account the potentially large, model-dependent corrections to m& and m T [35]). Moreover, 
a mass of m V3 ~ 0.05 eV is compatible with heavy Majorana masses ~ 10 14 GeV. Specifically, 
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Figure 1: The invariant Jcp as a function of the expansion parameter e for the case of SU(5) 
with large tan/3. The solid line was obtained numerically while the green dot-dashed (blue 
dashed) line uses only the first term (three first terms) in the power expansion in e. 



■ i i i i i i i i i i i i i i i i i i i i i i i i i i i i i ii I i i i i | i i i i | i i i i 




0* 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■* Q I 1 1 1 1 1 1 1 1 1 1 11 11 1 

0.5 1 1.5 2 2.5 3 1 2 3 

4> v „ (Rad) <^ v „ (Rad) 



Figure 2: The neutrino mixing angles (left) and the phases and values of the CP-violating 
Jarlskog invariant J C p (right), as functions of the phase 022 = ^ r 9( a 22)- The experimental 
limits on the three neutrino angles are denoted by horizontal lines. The values of the table 
match the values of the figure at Argfa^) = 0. 

for m Ua = 0.05 eV, the coefficients in Table [2] would predict M 3 = 2.62 x 10 14 GeV. We note 
the following points of interest. 
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(i) For almost all entries there is good agreement between the precise numerical values and the 
dominant term in the expansions. (In the Tables, for presentation purposes, we replace loge 
with its numerical value at e = 0.2 in the entries for leptogenesis and charged-lepton-flavour 
violation, which is a good approximation for the values of e of physical relevance). As seen in 
Fig. [T] (where we compare the numerical values of Jcp with those obtained by keeping the first 
and the three first terms in the expansions) the series expansion for Jcp is very accurate in the 
relevant range of e. 

(ii) The model prediction for #13 is numerically smaller than #12, but not parametrically smaller. 
Correspondingly, as seen in Fig. [2] (left), the model prediction for #13 lies not far below the 
present experimental upper limit. We also see in this plot that the model predictions for #23 
and #i2 lie generally near or within the bands allowed by experiment, whatever the value of the 
model phase parameter 022 ■ 

(iii) The model prediction for 8 is parametrically large and 0(1), as seen in Fig. [2] (right). We 
also see that the light-neutrino Majorana phases 0i t 2 are also generically large, and do not vary 
strongly with § v 22 - 

(iv) Observable CP-violation is expected in neutrino oscillations, since the Jarlskog invariant 
scales as e. Numerically, it is typically O(10~ 2 ), as also seen in Fig. [2] (right). 

(v) As seen in Table [31 leptogenesis is naturally embedded in the model, and is resonantly en- 
hanced via self-energy corrections, due to the quasi-degeneracy of the heavy Majorana neutri- 
nos: M2 ~ M3. In this example, the self-energy contributions are almost an order of magnitude 
larger than the vertex corrections. The resonant behaviour is also manifest in the expansions, 
which are divergent when keeping only the dominant terms (due to a term of the form 
which diverges for e —>■ —1/b). To demonstrate this, the expansions are displayed by writing 
separately the convergent (vertex) and the divergent (self-energy) contributions. As discussed 
above, in order o achieve successful leptogenesis, rfo must be in the range of eq. ljTTI) . For the 
values presented in Table [3], this implies a dilution factor of A ~ 500. Alternatively, by scaling 
m v D with global factor f s , we could allow for different values of the dilution factor A. In this 
case, the prediction m V3 ~ 0.05 eV could be preserved by also rescaling Mrr by a factor of 
f 2 . Then, eq. fflUj) indicates that e would be modified by a factor f 2 . For the set of textures 
discussed in this section, by setting f s = 0.0441 one could obtain r/ B in the range of eq. ffTTj) 
with A ~ 1 and M 3 ~ 5 x 10 11 GeV. In this case, the predictions for lepton flavor violation 
would be modified by a factor of approximately f 2 . For instance, the previous value of f s would 
imply a prediction of BR(/i — > ej) ~ 10~ 15 tanf3 2 which is still of the experimental interest for 
the large values of tan/3 assumed in this model, while the predictions for BR(r — > /ry) and 
BR(r — > cy) would become very small. 

(vi) As also seen in Table El the model predicts a large rate for BR(ii — > ej), close to the present 
experimental upper limit, as could have been expected for tan/3 > 40. This and the other rates 
for lepton-flavour violation have been estimated using (fl8l) . The values of the soft terms are 
such that the function J(mo, M\/ 2 , //) displays a small variation with m , as seen in Fig. 1 of 
[31] for A = 0, M 1/2 = 600 GeV, m = 300 GeV. In this case, T ~ 10- 18 G 2 F /a 3 . This function 
could increase by two orders of magnitude if My 2 were decreased, but this range of values is 
excluded by cosmology. Alternatively, the rate could be decreased by either increasing M1/2 
and mo, or by reducing the Dirac neutrino Yukawa couplings, with a corresponding adjustment 
of the heavy Majorana masses. 
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(vii) More unexpected is the suppression of BR{r — > /17), and especially of BR{r — > ej). These 
suppressions arise from the form of the charged-lepton mixing matrix V e which, in combination 
with the degeneracy of the heavy neutrinos, tends to align the rotated Yukawa matrix Y u in 
the base in which charged leptons are diagonal: Y u = V^Y v Un- We return to this issue in 
subsequent examples. 

3.3 Examples with Small tan/3 

In [32], the U(l) charges were chosen so as to cancel anomalies, and several solutions were 
found in the small-tan (3 regime of supersymmetric theories. Restricting their attention to 
solutions with C2 = C3 and textures where the heavy- Major ana mass matrix can be considered 
as diagonal, the authors of [32J found five different fits to SU(5) models, whose detailed formulae 
are given in Appendix I. The choices of U(l) charges correspond to different solutions of the 
anomaly-cancellation conditions, and the behaviours of the fits in (32] are very dependent on 
this. We distinguish: 

(a) Fits 1-3: In these fits, single right-handed neutrino dominance (SRHND) has been imposed 
on Mrr. This is achieved by choosing the charge n\ to be a negative number between —a/2 
and 0. Using the GUT values found in Table 11 of [32], we obtain the best fits when n\ is close 
to zero, as seen in Table HI 

(b) Fits 4,5: In these cases, we do not find fits with clear-cut SRHND. As one can see in TableHJ 
the required value for n\ is still zero 0. 

We have calculated the relevant observables, including the expected CP and charged-lepton- 
flavour violation, for representatives of these two classes of fits, in order to compare their 
behaviours for different sets of parameters within the same GUT framework. Since Fits 1-3 
and 4-5 have common characteristics, we focus on one solution from each group, and work 
specifically with Fits 2 and 4 of [32]. The values of e, o and n\ used in these fits are given in 
Table HJ and the chosen values of the numerical coefficients are tabulated in Tables [5] and 
[6] for Fits 2 and 4 of [32] , respectively. The corresponding predictions for neutrino masses and 
mixing angles, CP and lepton-flavour violation are given in Tables [7] and [H respectively. 

The values of tan j3 for Fits 2 and 4 are obtained by using the Yukawa couplings Y e and the 
experimental value of m T ; these lead to tan/3 = 2.04 for Fit 2 and tan/3 = 5.15 for Fit 4. The 
textures Y v and M RR are scaled so that m U3 = 0.05 eV when the largest right-handed neutrino 
mass becomes M 3 = 5 • 10 14 GeV. For Fit 2, this implies that 

Y U ->Y V - 0.46 • e- 9/2 , Mrr -> Mrr/ (E) • 2.1 • 10 14 • e ~ 19/2 , (31) 

while for Fit 4 

Y V ^Y V - 0.42, M RR -> Mrr/ (E) ■ 2.34 • 10 14 . (32) 
The expansions of the observables in powers of e are obtained as discussed above. 



In principle, one could also investigate what happens for other forms of Man, since in general there might 
be dominant off-diagonal terms (depending on the singlet charges). However, such a study would go beyond 
the scope of this paper. 
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Figure 3: The Jarlskog invariant Jcp as a function of the expansion parameter e for two SU(5) 
models at low tan/3. The textures correspond to Fits 2 and 4 of [32] for the values of the 
coefficients in Tables^ and\£h The solid line was obtained numerically while the dot-dash 
(dash) line uses only the first term (three first terms) in the power expansion in e. 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 



4> V „ (Rad) c^ v „ (Rad) 



Figure 4: The neutrino mixing angles $23,12,13 i n various SU(5) models, as functions of the 
phase = Arg^a^- The textures correspond to Fits 2 and 4 of fSBj for the values of the 
coefficients in tables and [3 The experimental limits on the three neutrino angles are 
denoted by horizontal lines. 

As in the first model example, the changes induced by altering the numerical values of the 
coefficients affect the coefficients in the expansions of the mass matrices, but not the powers 
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4> v „ (Rad) (j) v „ (Rad) 



Figure 5: The neutrino phases and values of the CP-violating Jarlskog invariant Jcp as func- 
tions of the phase (p 22 = Arg(a 22 ) for Fit 2 and Fit 4, respectively. 



Case 


e 


a 


m 


Fit 2 


0.2 


19/2 





Fit 4 


0.19 


0.5 






Table 4: Parameter values obtained for Fits 2 and 4 of [32], corresponding to the textures of 
eqns. (9.2, 9.3, 9.4) of that paper. The values are specified at the GUT scale, the stand for 
the flavour charges of the right-handed neutrinos and o is the flavour charge of the singlet field 
that, through its vev, generates the right-handed neutrino mass hierarchies. 





Fit 2 of [32J 


Charged leptons 


a\ 2 = 0.55, a e 21 - -0.6, of 2 = 2, a e 23 = -0.7e j7r / 4 , a\ 2 = 1.5, o| 3 = 1.4 






M RR 


a" = 1.1, of 2 = 1.2, of 3 = 1.2, a l 2 \ = 2, ag = 0.45 



Table 5: Choice of coefficients to reproduce the fermion observables for Fit 2 of J3J% . Coefficients 
not appearing in the Table are set to unity. 





Fit 4 of [32J 


Charged leptons 


af 3 = 0.8, a\ 2 = -2, a\ 3 = 0.8e^ 5 , a e 33 = 1.1, 


m u D 


a\ 3 = 0.5, a» 22 = 0.7e-/ 4 , a\ 2 = 1.9, a? m = 0.7 


M RR 


£^ 2 = 2,a^ = 0.65,a^ = 1.4 



Table 6: Choice of coefficients to reproduce the fermion data for Fit 4 of IWj . Coefficients not 
appearing in the table are set to unity. 
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Observables 


Series Expansions 


Numerical Values 


m U3 (eV) 


-0.055 + 0.0071e- 3 / 2 + 0.0084e~ 5 / 4 


0.05(-0.055) 


m Ul /m Us 


2.21e 2 + 0.i2e 5 / 2 -0.051e n / 4 


0.16(0.09) 


m U2 /m U3 


4.13e 3 / 2 - 0.29e 2 - 13.82e 9 / 4 


0.20(0.37) 


M x /M 3 


0.82e - 0.98e 5 / 4 + 0.55e 3/2 


0.06(0.16) 


M2/M3 


0.636 1 / 4 - 1.49eV 2 + 1.73e 3 / 4 


0.12(0.42) 


m Uee (eV) 


0.017e 4 / 3 + 0.04e 5 / 3 - 0.85e 2 


0.0042(0.002) 


6*23 


1.08 + 0.17e 3 / 4 -0.012e 4 / 3 


0.91(1.13) 




0.153c 1 / 4 - O.lOe 2 / 3 + 0.066e 3 / 4 


0.56(0.10) 


#13 


0.21e 2 / 3 + 0.25e + 0.0004e 4/3 


0.21(0.07) 


s 


-.45 + 0.059e 1/3 - 0.31e 5/12 


0.43(-0.45) 


Jcp 


0.0029e 1 + 0.003e s/ + 0.0009c' 7 


0.0088(0.00066) 


J. 

0i 


010 i r\oo 1/3 c\ 5/12 

2.13 + 0.23e ' — O.ooe ' 


2.94 (2.13) 


02 


2.04 + 0.23e 1 / 3 - 1.25e 5 / 12 


4.6 (2.03) 


Cl 


0.044e 3 / 2 - 0.079e 7 / 4 - O.OOle 9 / 4 ) 


0.00046(0.0039) 


62 


0.00046e 3 / 4 + 0.0085e + 0.07e 5 / 4 


0.0007(0.00014) 


€3 


0.00005e - 0.06e 5 / 4 + 0.2e 3 / 2 


0.0003(0.000001) 


Vb x A 




1.18 x 10~ 7 (3.11 x 10- 7 ) 




0.11e 2 / 3 + 0.023e n / 12 -2.64e 


0.44(0.037) 






8.5 x 10~ 12 (6.1 x 10" 14 ) 


x m| 


0.057e 2 / 3 + 0.012e n / 12 - 0.58e 


0.11(0.02) 






5.95 x HT 13 (1.74 x 10~ 14 ) 


x mf 

re L?-} 


0.24 + 0.05c 1 / 4 + 0.002c 1 / 2 


0.34 (0.24) 


BR(r -> /i 7 ) 




5.1 x 10~ 12 (2.6 x 10~ 12 ) 



Table 7: Values of observables predicted by Fit 2 of f32^ . for tan/? = 2.04. The rate of con- 
vergence of each expansion in e can be judged from the relative magnitudes of the expansion 
coefficients and by comparing the exact numerical values of the observables with those obtained 
by keeping only the first term in each expansion (in parenthesis). 
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Observables 


Series Expansions 


Numerical Values 


m V3 (eV) 


-0.054e + 0.024c- 1 / 2 + 0.02e~ 3 / 8 


0.05(-0.054) 


m vl /m V3 


0.26e 3 / 2 + 0.50e 7 / 4 + 0.34e 15 / 8 


0.069(0.022) 




1.046 1 / 2 - 1.94e 3 / 4 + 1.24e 7 / 8 


0.21(0.45) 


M x /M 3 


0.046 - 0.046 1 / 8 - O.OOSe 1 / 4 


0.034(0.046) 


M 2 /M 3 


0.95 - 1.64c 1 / 8 + 1.41c 1 / 4 


0.15(0.95) 


rn Vee (eV) 


0.026e 7 / 4 + 0.15e 2 + 0.05e 17 / 8 


0.0027(0.0014) 


#23 


0.90 - 0.027e 3/8 - 0.53e 1/2 


0.87(0.66) 


6*12 


0.56 1 / 2 + 0.94e 3 / 4 + 0.026e 7 / 8 


0.53(0.22) 


6*13 


0.54e - 0.27e n / 8 - 0.88e 3 / 2 


0.04(0.10) 


5 


2.56 - 0.025e 1/4 + 0.03e 3/s 


4.87(2.56) 


J 

Jcp 


0.036e J/ + 0.068e (/ — 0.017e 7 


r\ r\r\ a { r\ r\r\^>\ 

0.004(0.003) 


J. 

01 


cr in 1/4 i nnn/^ 3/8 

5.19 — 0.074c' + 0.096e 1 


6.17(5.20) 


02 


2.056 + 0.0256 1 / 4 + 0.096e 3 / 8 


2.91(2.06) 


ei 


0.000085e 3 / 8 + 0.00017c 1 / 2 - 0.00013e 5 / 8 


0.00087(0.000046) 


C2 


0.0013e 3/8 - 0.00126 1 / 2 + 0.00053e 5/8 


0.0012(0.00071) 


£3 


0.0014e 3 / 8 + 0.00136 1 / 2 + 0.0006e 5 / 8 


0.0002(0.00074) 


77b x A 




1.62 x 10~ 7 (1.03 x 10" 7 ) 


16^ x m 2 


0.17e 5 / 8 + 0.00024e 3 / 4 - 0.046e 7/8 


0.05(0.06) 






7.8 x 10~ 13 (1.04 x 10~ 12 ) 


l&r" x m 2 
« Lis 


0.32e 5/8 + 0.0066e 3/4 - 0.062e 7/8 


0.16(0.12) 


RR(r -> e7) 




7.9 x 10~ 12 (3.9 x 10" 12 ) 


x mf 


0.24 + 0.126 1 / 8 + 0.0356 1 / 4 


0.48(0.24) 


BR(r -> // 7 ) 




6.6 x 10 -u (1.7 x 10~ n ) 



Table 8: Values of observables predicted by Fit 4 of %3E/, for tan/5 = 5.15. The level of 
convergence of each expansion in e can be judged from the relative magnitudes of the expansion 
coefficients and by comparing the exact numerical values of the observables with that obtained 
by keeping only the three first terms in each expansion (shown in parenthesis). 
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of e. However, when varying the absolute values and the phases of the coefficients afj, the 
numerical values of the observables are generally affected. We note the following points in 
connection with these two examples. 

(i) In both fits, the dominant terms in the expansions generally agree only in order of magnitude 
with the exact numerical values, for the values of e required to reproduce the correct fermion 
mass hierarchies. The leading term in many expansions is of order z = e 1 / 8 , and yields good 
numerical results only when z ~ 0.5 — 0.75 (i.e., e ~ 0.05 — 0.1). This slow convergence leads 
to the numerical instabilities that are evident in many observables, for example in Jcp, as 
shown in Fig. [31 These instabilities may be introduced by combinations of the following: (a) 
the powers of e in the series are often very close to each other (e.g., ae 3//8 + 6e 4//§ + ce 5//8 ), and 
(b) the relative signs of the terms in the expansions may change with the choice of a^. 

(ii) The model predictions for 6*23 is 0(1), and that for #12 is suppressed only by C(e 1//4 ). As seen 
in Fig. HI the numerical values of these angles are reasonable for 0^ in Fit 2, in which case the 
prediction for #13 is close to the present experimental upper limit, even though parametrically it 
is suppressed by 0(e 2 / 3 ). In the case of Fit 4, the values of #23,12 are reasonable for 0.3 < 022 < 1> 
but #13 is always considerably below the present experimental upper limit. 

(iii) The complex phases 5 and 0^2 are 0(1) in both fits, but vary widely with 022' as seen i n 
Fig. E 

(iv) Observable CP-violation is to be expected in neutrino oscillations, although the Jarlskog 
invariant Jcp is significantly smaller than in the case of large tan/5, as also seen in Fig. [51 
particularly for Fit 4. In this later case, there is a rather good agreement between the numerical 
value and the dominant terms in the expansion. 

(v) Both Fit 2 and Fit 4 may accomodate successful leptogenesis, and also lead to lepton flavour 
violation that is close to the present experimental upper limits. Values of rjb in the range of 
eq. (|17[) would be obtained with the values presented in Table [7] and Table [S] if the dilution 
factor A ~ 200. As we discussed in the previous subsection, a global scaling factor f s for m v D 
and the corresponding rescaling of Mrr would imply a modification of the leptogenesis (and 
lepton-flavour violation) predictions. We find that a factor f s ~ 0.06 — 0.07 would predict tjb in 
the range of eq. (fTTj) with A ~ 1 and M 3 ~ 10 12 GeV. However, this would imply a reduction 
of about two orders of magnitude in the predictions for charged-lepton-flavor violation. 

(vi) The numerical results shown in Fig. H] further demonstrate the inter-correlations between 
the different physical parameters, indicating how the experimental limits on neutrino masses 
constrain the allowed range for the angles, phases, and Jcp- We see that, even for the same 
GUT group and range of tan j3, the relations between the observables display certain differences, 
reflecting among others the role of the coefficients in obtaining viable solutions in these schemes. 

3.4 Flipped SU(5) Model 

In the case of the flipped SU(5) GUT model, the fields Qi,d^ and v\ of each family belong 
to a 10 representation of SU(5), the u\ and Lj belong to 5 representations, and the e\ fields 
belong to singlet representations of SU(5). Thus the correlations between the U(l) charges 
of the different matter fields are different from those in conventional SU(5). These particle 
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assignments imply a symmetric down-quark mass matrix, whereas the structure of the up- 
quark mass matrix depends on the charges of the right-handed quarks. However, as these are 
the same as the charges of the left-handed leptons, the up-quark mass matrix is constrained by 
the need to generate large mixing for atmospheric neutrinos. 

The simplest example that matches the charged-fermion mass hierarchies (adjusting coefficients 
so as to match the experimental value of V^) is [16] 



Once again, the form of the heavy Majorana mass matix depends on the charge of the field 
S. For instance, for a = 0, it will be similar in structure to the down-quark mass matrix. 
The contribution from the up-quark sector to is generically small in this model, leading to 
V c b — \Jm s /m h . (indeed, for the up-type quark hierarchies, it turns out that e ~ 0.2, e ~ e 4 
and \y\ ~ 2). This is too large, and requires a significant adjustment of the 0(1) coefficients. 

However, the problem with the large value of V c b can be avoided by combining flipped SU(5) with 
a non-Abelian flavour group, or by adding a second singlet field with different transformation 
properties under the flavour group. In that way, one could obtain solutions similar to those 
of [IB|, with 



However, even after overcoming this obstacle, it is very difficult to obtain naturally viable 
solutions: if e ~ e 4 , the Dirac-mass hierarchies turn out to be too large to generate interesting 
solutions. Of course, the fact that the see-saw mechanism is hard to implement does not mean 
that an appropriate m e ff could not be generated by alternative mechanisms. However, the 
symmetries require that 



indicating that, unlike in conventional SU(5), it is difficult to generate large solar mixing. 

As in the case of conventional SU(5), one could either (i) study the most generic cases with 
half-integer charges, requiring the (3,3) flavour charges to be non-zero, and using the see-saw 
conditions to obtain large mixing for the solar neutrinos, or (ii) add further fields and flavour 
groups. However, any departure from the simple structure described so far introduces significant 
model dependence, and results in a loss of predictivity, unless this additional structure is 
predicted (and constrained) by concrete theoretical considerations. In the subsection that 
follows, we discuss what may be the most natural way to proceed, if one believes in the existence 
of an underlying string theory. 




(33) 




(34) 




(35) 
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3.5 String-Derived Flipped SU(5) Models 



So far, we have been discussing models based on a single U(l) flavour symmetry, and with only 
one field 9 whose vev provides a single expansion parameter in the mass matrices. However, 
in realistic models this need not be the case. On the contrary, one may expect several fields 
to be involved in the generation of mass terms, while additional flavour symmetries may be 
relevant. While this seems to limit predictivity, in models based on an underlying string theory, 
the string symmetries (translated into selection rules) impose strong constraints on the mass 
and mixing matrices. 

In the previous sections we observed that there is significant model-dependence in the results, 
and that the values of certain coefficients must be in specific ranges in order to give viable 
solutions. Optimally, we would like to understand these values from prior principles. In a 
realistic model, they could be related to the vacuum expectation values of fields that are 
constrained, for example, by considerations on flat directions. 

In specific string-derived models, due to the string selection rules, several texture zeros are 
generated in the matrices, and these could in principle lead to strong constraints on certain 
observables. Indeed, the constraints are so strong that one could expect that it might be 
difficult to generate the required entries that fully reproduce the observed fermion patterns. 
Nevertheless, it was found in [361 [37] that it is possible to accommodate all data in a natural 
and generic way within a flipped SU(5) x U(l) string model [38]. Relevant aspects of this model 
are reviewed in Appendix II: the theory contains many singlet fields, and the mass matrices 
depend on the subset of them that get non-zero vev's, i.e., on the choice of flat directions in 
the effective potential. 

The flat directions and the quark masses and mixings for this model have been studied in [36J, 
where it was found that 
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A 2 A 3 <5 23 
V A 5 A 3 3 

As already remarked, the relevant field definitions are given in Appendix II. For instance, A2A5 
is a combination of hidden-sector fields that transform as sextets under SO (6). 

The lepton mass matrices were studied in [37], where it was shown that: 
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where 



y = F 5 ,r = 03, 6 = $ 3 i$ 3 i0 2 , c = $23, 03 = 



(40) 



Note that these matrices are given in the flipped SU(5) field basis. However it is easy to pass to 
the weak interaction eigenstates by an appropriate rotation. For instance, the weak-interaction 
eigenstates for the light neutrinos have the following assignments: 

ve - h - 0(A 2 A 5 )/ 2 , u„ - h + 0(A 2 A 5 )/ 5 v T -> A, (41) 

while the flipped SU(5) basis (/s, / 2 , fi) is related to (J4"lj) by the rotation 



l/ m 



/ 1 - i(A 2 A 5 ) 2 A 2 A 5 

-A 2 A 5 1 - |(A 2 A 5 ) 2 I . (42) 

V 01 



The forms of the mass matrices depend on the various field vev's. The analysis of quark masses 
pointed towards x = A 2 A 5 = 0(1) (and potentially large solar mixing, already from the 
charged-lepton sector), as well as a rather suppressed value of / = 4 <C 1 (since V c b ~ A 2 A 5 4 ). 
Moreover, from the analysis of flat directions [35], it was concluded that $ 3 i$ 23 = 0(1) is large. 
In addition, the flatness conditions [3B] relate $ 3 i, $ 3 i and 2 , and can be satisfied even if all 
the vev's are large, as long as $ 3 i<I> 3 i and $ 23 $ 23 are not very close to unity. As for the 
decuplets that break the gauge group down to the Standard Model, we know that the vev's 
should be ~ Mgut/M s . In weakly-coupled string constructions, this ratio is ~ 0.01. However, 
the strong-coupling limit of M theory offers the possibility that the GUT and the string scales 
can coincide, in which case the vev's could be of order unity. 

The 2x2 form of the charged-lepton mass matrix in fact puts severe constraints on the fields 
involved. On the other hand, flatness conditions and quark masses do not give any information 
on the vev of the product T 2 T 5 . However, this combination is to some extent constrained by 
the requirements for the light neutrino masses [37]. Finally, the field 3 is the one for which 
we seem to know least and we can consider two extreme possibilities: a very large value 0(1) 
and a very small value. For large 3 , if 4 ~ F 5 , F\ as would be expected in weak-coupling 
unification schemes, the entries of Jn e // are & 11 of the same order of magnitude, and large — v e 
and Up — v T mixings are naturally generated via the neutrino mixing matrix. 

Charged-lepton mass hierarchies pose severe constraints on the model parameters. From ([3 
we get 

f 
u' 2 (l + x 

while / and x are also constrained from the quark masses, by 



m c /m t ~ / x 0(1), (44) 
V cb ~ xf ~ 0.044. (45) 



Therefore, if / is 0(0.01), x can be 0(1) while the charged leptons mass ratio can be satisfied 
with u of 0(0.1). 
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The charged-lepton mixing matrix V e is obtained by the diagonalization of mg in the weak- 
interaction basis: 

Vl{Vr L ) T mm\{VrjV: = Diag{mlm%ml). (46) 
We parametrize the product of these two rotations by defining: 

/ l/N g/N \ 
Qe = V e Vr L ~ -9/N l/N , (47) 

V o o i) 

where g is a mixing parameter to be determined by the neutrino mass matrix, and N is a 
normalization factor for the sector 1-2: N = 1/(1 + |g| 2 ) 1//2 . 




Figure 6: The neutrino mixing angles (left) and the phases and value of the CP-violating 
Jarlskog invariant Jqp (right), as functions of the absolute value of g. The experimental limits 
on the three neutrino angles are denoted by horizontal lines. 



Passing finally to the neutrino mass matrix, we see the following. In the basis where the 
charged-lepton masses are diagonal, (i.e. — > Q T e ■ m„), the vev's for the remaining fields 
and M VR enter in m K ff as combinations: 
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eff 



/ (Gn + p-HnVN 2 (G 12 + p ■ H 12 )/N 2 
(da + p • H 12 )/N 2 (G 22 + p • H 22 )/N 2 
V q ■ (1 + gx - 2x 2 )/N -q • (x - g(l - 2x 2 ))/N 



q- (1 + gx- 2x 2 )/N 
-q-(x-g(l-2x 2 ))/N | . (4? 
4q 2 • x 2 



Here, p, q are combinations of vev's defined as: 



P 



c • r ■ t y • f i • r 

while the functions Gij(x,g), Hij(x,g), are of 0(1): 
G n (x,g) = g 2 -2gx + x 2 , 



(49) 



23 



G 12 (x,g) = x - g 2 x + g(-l + x 2 ), 

G 22 (x,g) = l + 2gx + g 2 x 2 , 

H n (x,g) = 3g 2 x 3 + x(-l + Ax 2 ) -2gx(x + 2x 3 ), 

Hi 2 {x,g) = gx(-l + x 2 ) + x(x + 2x 3 ) - g 2 x(x + 2x 3 ), 

H 22 {x,g) = 3x 3 + g 2 x{-l + Ax 2 ) + 2gx{x + 2x 3 ). (50) 



The phase of q can be eliminated with a rotation g 3 ■ m e ff ■ q^, where q$ = diag(l, 1, Arg(q)/2). 
The neutrino data can be fit by choosing suitably the parameters p, q, g and x. The results 
presented in Table M correspond to following choice of parameters: 

x = 0.8 • e a3 \ g = 0.7 • e" a94 \ p = 0.43 • e 2 ' 2i , q = 0.25, (51) 

while in Fig. Owe vary the value of \g\ while keeping constant the rest of the parameters. 

To estimate the values of the observables that involve the right-handed neutrinos, we need 
the individual values of the vev's contained in p, q. Since, b, c, r are 0(1), t ~ f 3 and 
also y ~ Fi ~ /. In order to provide some numerical values we can make the assignment 
y = Fi = t 1 / 3 = |/| = 0.01. Then, from the values of p, q we can get: 

\r/b\ = 1/4, \c/b\ = 1.84. (52) 

Despite the phase dependence of M v r being rather complex since all the vev's contain phases, 
these appear on the observables e\^ 2 and BR(/i — > e^y) in the same combination as the phase 
of p. For field vev's chosen as above, the matrix M VR in (|39l) contains two almost degenerate 

eigenvalues M 12 ~ \{M VR )\ 2 + (Mur)^ 1 ^ ' ■ The predictions of the string-derived flipped SU(5) 
model for the neutrino masses, mixing angles, CP and charged-lepton flavour violation are 
shown in Table [9j 

The following are some notable features. 

(i) Since the vev's are fixed in the string-derived flipped SU(5) model, there is no auxiliary 
expansion parameter, as in the previous SU(5) models. However, it is interesting to study 
observables function of g, as shown in Fig. El 

(ii) We see in Fig. [6] that #23 is large and within the experimental range for all values of g, 
whereas 8\ 2 lies within the allowed range for \g\ > 0.5, and #13 is acceptably small for \g\ > 0.6 
and may lie close to the experimental upper limit. 

(iii) We also see in Fig. [6] that 5 and the two low-energy Majorana phases are expected to be 
quite large. 

(iv) We also see that Jcp may be quite large, though there is also the possibility that it might 
(be close to) vanish(ing). 

(v) We see in Table M that leptogenesis is quite possible. Values of r\ h in the range of eq. ffTTI) 
with the values presented in Table [9] are consistent with a dilution factor A ~ 80. In this case, 
we find that a global scaling factor f s ~ 0.11 on m v D would predict tjb in the range of eq. (jTTj) 
with A ~ 1 and M 3 ~ 10 12 GeV. Such a scaling would decrease the prediction for BR(/i — > ej) 
by two orders of magnitude, making it compatible with the experimental bound for a lower 
range of supersymmetric masses. 
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Observables 


Numerical Values 


m V3 (eV) 


0.056 


m Ul /m V3 


0.008 


m u ,Jm U3 


0.17 


Mi/Mz 


1.4 x 10" 4 


M 2 /M 3 


0.98 


m Uee (eV) 


0.0027 


#23 


0.84 


012 


0.61 


&13 


0.15 


5 


4.1 


Jcp 


0.014 


01 


2.86 


02 


6.26 


ei 


1.5 x 10~ 7 


£2 


3.1 x 10" 3 




3.2 x 10" 3 




4.5 x 10" 8 


BR(fi -> e 7 ) x (ta ^ )2 


2.4 x 10~ 1U 



Table 9: Values of observables predicted in the string-derived flipped SU(5) model for the choice 
of parameters given in the text. 



(vi) The values of Table [9] correspond to M3 = 10 14 GeV. However, these imply a BR(/i — > ej) 
above the experimental bound (unless the supersymmetry breaking masses become very large), 
due to 0(1) entries in the Dirac neutrino mass matrix (1351) . Clearly, by lowering the heavy 
right-handed neutrino masses, and thus also the Dirac neutrino couplings, the rates are altered 
accordingly. For example, M 3 = 5 x 10 12 GeV would yield BR(/i — > cy) ~ 1.4 • 10~ 12 • tan/3 2 , 
while £1 ~ t2 ~ 1.5 • 10 -4 . The expressions for ei^ and BR(/i — > ej) contain Y v in a basis where 
M UR is diagonal, so their values have the same phase dependence as in the 1-2 sector of m e jf. 
Hence, with the choice of parameters of (IBTj) . the results presented in Table [9] are independent 
of the individual phases of the parameters included in the definition of p in (1491) . 

(vii) On the contrary, due to the lack of mixing in the third generation in this model, r decays 
are not observable in this case. 



4 Conclusions 

We have studied in this paper the predictions for CP and charged-lepton-number violation 
in different SU(5) models, including a string version of flipped SU(5). Because of the inter- 
correlations between the different neutrino observables, it was in each case possible to obtain 
quite specific predictions, which differ significantly between the models. 



25 



The general indications are that accessible rates for CP and charged-lepton-number violation are 
to be expected. However, due to the different expectations for the couplings and angles in the 
various models, the magnitudes of the observables may be quite different, even between models 
based on the same group. We note that #13 may be close to the present experimental upper limit 
in some models, but much smaller in others. Also, the predictions for the Jarlskog invariant 
J cp vary significantly. Among the most sensitive observables are those for charged- lepton- 
flavour violation and leptogenesis, and cases with a resonant enhancement of the expected 
lepton asymmetry for degenerate heavy Majorana neutrinos are particularly interesting. 

It is crucial to keep in mind that the various models sometimes contain instabilities due to 
cancellations, to which the results may be sensitive. In principle, for cases such as the flipped 
SU(5) string model, we expect that the coefficients are to a good extent predicted by the vev's 
of the fields. However, one should not forget that there is a sum of non-renormalisable terms, 
implying that, if the higher-order terms involve fields with high vev's, they can effectively 
introduce multiplicative factors, particularly in the case of many fields. In principle, once 
some of the observables discussed here are measured, one may use the correlated data to refine 
the predictions with better knowledge of the possible ranges of coefficients and the level of 
tuning via sub-determinant cancellations. However, such an analysis is premature ahead of the 
corresponding measurements, and would go beyond the scope of this paper. 



Appendix I: GUT Fits 2 and 4 of p2| 

We give here the matrix parametrizations found in Fits 2 and 4, as used in the text. 
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Fit 4: 
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Appendix II: Flipped SU(5) Particle Assignments 



In this Appendix we tabulate for completeness the field assignment of the 'realistic' flipped 
SU(5) string model [38], as well as the basic conditions used in [36] to obtain consistent flatness 
conditions and acceptable Higgs masses. 



^(10, 0,0,0) 
F 2 (10,i, 0,-1,0,0) 
F 3 (10 ,i 0,0,|,-|) 
F 4 (10,|,-|, 0,0,0) 
F 5 (T0,-i,0,i,0,0) 


A(5,-|,-i, 0,0,0) 
/ 2 (5,-|,0,-i, 0,0) 
/3(5, —5, 0, 0, i, \) 
/ 4 (5, |, |,0,0,0) 
/ 5 (5, -|,0, -|,0,0) 


3(1, |, -±, 0,0,0) 

3(1, §,0,-1, 0,0) 

l 3\ L ^ 2' U ' U ' 2' 2> 

£1(1, -|, |,0, 0,0) 
3(l,f,0,-±,0,0) 




^(5,-1,1,0,0,0) ^(5,-1,0,1,0,0) 
^ 5 (5,-l,-i-i0,0) 


^(5,-1,0,0,1,0) 



045(1,0, 5, 


5,1,0) 


+ (l, 0,1, -1,0,1) 


0.(1,0,1,-1,0,-1) 


$23(1,0, 0, 


-1,1,0) 


$31(1,0,1,0,-1,0) 


$i 2 (l, 0,-1, 1,0,0) 


0(1,0,1,- 


5,0,0) 


$(1,0,0,0,0,0) 





A^O, 1,6, 0,-1,1,0) 
A 4 (0, 1,6, 0,-1, 1,0) 



A 2 (0,l,6,-i,0,l,0) 
A 5 (0,l,6,l,0,-1,0) 



A 3 (0, 1,6, -1,-1, 0,1) 



T 2 (0, 10, 1,-|,0, 1,0) 
T 5 (0, 10, 1,-1, 0,1,0) 



T 3 (0, 10, 1,-1,-1, 0,1) 



71(0,10,1,0,-1,1,0) 
T 4 (0, 10, 1,0, 1,-1,0) 



-, 2' 2' ' 2' 



Table 11: The chiral superfields are listed with their quantum numbers 138]. The Fi, fi, i\, as well 
as the hi, hij fields and the singlets are listed with their SU(5) x U(l)' x U(l) quantum numbers. 
Conjugate fields have opposite U(l)' x U(l) 4 quantum numbers. The fields Aj and Ti are tabulated in 
terms of their U(l)' x SO(10) x 50(6) x U(l) 4 quantum numbers. 

As can be seen, the matter and Higgs fields in this string model carry additional charges under 
additional U(l) symmetries [38]. There exist various singlet fields, and hidden-sector matter 
fields which transform non-trivially under the SU(4) x SO(10) gauge symmetry, some as sextets 
under SU(4), namely Ai !2 ,3, 4 ,5, and some as decuplets under SO(10), namely Ti i2 ,3, 4 ,5. There 
are also quadruplets of the hidden SU(4) symmetry which possess fractional charges. However, 
these are confined and do not concern us further. 

The usual flavour assignments of the light Standard Model particles in this model are as follows: 

fi:u,r, f 2 :c,e/fx, / 5 :t,/i/e, 
F 2 : Q 2 , s, F 3 : Q 1 , d, F 4 : Q 3 , b, 



: r. 



-2 ■ e , 



/'■ 



(55) 



27 



up to mixing effects which are discussed in more detail in |36j. We chose non-zero vacuum 
expectation values for the following singlet and hidden-sector fields: 

$31, $31, $23, $23, 02, 03,4, <t>~ , + , 045, 045, A 2 , 3 ,5, ^2,4,5- (56) 

The vacuum expectation values of the hidden-sector fields must satisfy the additional con- 
straints 

T 3 2 4j5 = Tj • T 4 = 0, A| 5 = 0, T| + Aj = 0. (57) 
For further discussion, see |36j and references therein. 
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